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Quantum phases of Bose-Einstein condensates with synthetic 
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Chunlei Qu, Kuei Sun, and Chuanwei Zhang* 

Department of Physics, The University of Texas at Dallas, Richardson, Texas 75080-3021, USA 

The experimental realization of emergent spin-orbit coupling through laser-induced Raman tran¬ 
sitions in ultracold atoms paves a way for exploring novel superfluid physics and simulating exotic 
many-body phenomena. A recent proposal with the use of Laguerre-Gaussian lasers enables another 
fundamental type of coupling between spin and orbital angular momentum (SOAM) in ultracold 
atoms. We hereby study quantum phases of a realistic Bose-Einstein condensate (BEC) with this 
synthetic SOAM coupling in a disk-shaped geometry, respecting radial inhomogeneity of the Raman 
coupling. We find that the experimental system naturally resides in a strongly interacting regime 
in which the phase diagram significantly deviates from the single-particle picture. The interplay 
between SOAM coupling and interaction leads to rich structures in spin-resolved position and mo¬ 
mentum distributions, including a stripe phase and various types of immiscible states. Our results 
would provide a guide for an experimental investigation of SOAM-coupled BECs. 
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I. INTRODUCTION 

The interplay between a single particle’s spin and or¬ 
bital motion, or spin-orbit coupling, plays a crucial role in 
various nontrivial many-body phenomena, such as topo¬ 
logical insulators and superconductors [1, 2] in fermionic 
systems, as well as exotic spinor condensates and super¬ 
fluids in bosonic ones. Considerable effort [3-40] has 
been devoted to the investigation of spin-orbit coupled 
many-body physics in ultracold atoms due to the sys¬ 
tem’s high tunability, disorder-free environment, and, 
most importantly, synthetically inducible spin-orbit cou¬ 
pling through the atom-light interaction [41-43]. The 
recent experimental realization [3-15] of spin and lin¬ 
ear momentum (SLM) coupling with the use of counter- 
propagating Gaussian lasers to induce two-photon Ra¬ 
man transitions [41] enables the direct observation and 
manipulation of SLM coupled Bose-Einstein condensates 
(BECs) and degenerate Fermi gases. In spin-orbit cou¬ 
pled BECs, interaction effects are expected to be es¬ 
sential for newly emergent quantum phases [16]. How¬ 
ever, compared with the typical kinetic energy scale in 
the current experimental system, i.e. the recoil energy 
El = h 2 /2M\ 2 with atomic mass M and the laser’s 
wavelength A ~ 1/im, the interaction effects are relatively 
weak such that most behaviors of the SLM coupled BECs 
show little difference from the single-particle ones [3]. 

The realization of another fundamental type of spin- 
orbit coupling, namely the spin-orbital angular mo¬ 
mentum (SOAM) coupling, in ultracold atoms has 
been proposed in our recent work [44], with the use 
of co-propagating higher-order Laguerre-Gaussian (LG) 
lasers [44-46] (see Fig. 1). The two LG laser beams 
with different orbital angular momentum (OAM) stim¬ 
ulate the Raman transitions between two atomic hyper- 


* Author to whom all correspondence should be addressed: chuan¬ 
wei. zhang@utdallas.edu 


fine states, \f) and |^), and concurrently impart spin- 
dependent OAM to the atoms [47-49]. In such a system, 
the rotational motion defines the kinetic energy scale 
Eq = h 2 /2MR 2 with R ~ 15/im of the order of the LG 
laser’s beam waist, which is much smaller than the recoil 
energy El and thus greatly enhances the many-body ef¬ 
fect. Therefore, the SOAM coupled system provides an 
opportunity to explore the spin-orbit coupled physics in 
a strongly interacting limit. A SOAM coupled ring BEC 
has been analyzed to reveal several salient features of 
the aforementioned physics [44]. However, the effect of 
the LG laser beams’ radial distributions is integrated out 
and plays no role in the ring system. This inhomogeneous 
beam intensity will generate radial-dependent potential 
as well as radial-dependent SOAM coupling in a two- 
dimensional (2D) system, and it will result in an inter¬ 
esting energetic competition against kinetic energy and 
interaction. 

In this paper, we study the quantum phases of a realis¬ 
tic 87 Rb BEC in a disk-shaped geometry with SOAM cou¬ 
pling. We first discuss the single-particle physics of the 
2D system and compute interacting phase diagrams as 
a function of Raman coupling and detuning with experi¬ 
mental parameters. We then analyze the detailed wave- 
functions of each individual phase and show their rich 
structures attributed to the 2D geometry and interaction. 

(a) (b) 



FIG. 1. (Color online) (a) Illustration for generating SOAM 
coupling and (b) Raman transitions induced by the LG laser 
beams. 
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In Sec. II, we build the model Hamiltonian for a SOAM 
coupled BEC. The single-particle physics is discussed in 
Sec. III. Section IV presents the main results. We also 
present the phase diagram of the interacting BEC and 
discuss the individual phases. In Sec. V, we extend the 
investigation to the trapping potential and LG beams’ 
OAM-induced quantum phase transitions. Section VI 
contains the conclusion. 


II. MODEL HAMILTONIAN 

We consider a pair of co-propagating LG beams of the 
semiclassical form (j = 1 , 2 ) 


r) = 
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shining on a disk-shaped BEC, where (r, </>, z) are 
the three cylindrical coordinate variables as shown in 
Fig. 1 (a). For simplicity, we assume that the two LG 
beams have the same amplitudes f^o and waist re, but 
possess opposite OAM li = +Z and I 2 = —l- We regard 
the two atomic hyperfine states as two pseudo-spins ||) 
and ||), respectively. Similar to the experiment for SLM 
coupling [3], an external Zeeman field B induces a split¬ 
ting huz between the two states, which also defines the 
frequency difference of the two LG lasers Scu. The two 
pseudo-spins are then coupled through a two-photon Ra¬ 
man transition as illustrated in Fig. lb. Both LG beams 
are assumed to be red-detuned from the excited state by 
A < 0 , while S = Scu — hwz is the two-photon Raman 
transition detuning. 

The equilibrium and dynamic properties of the BEC 
with LG beams induced SOAM coupling are governed 
by the Gross-Pitaevskii (GP) equation 
where T = (7/q,Vq) T is the spinor wave function of the 
BEC. The GP Hamiltonian 


HgP= 2M + ^ T(r) + ^ LG(r) + ^ l(r) 


( 2 ) 


comprises three potential parts: ( 1 ) the external trapping 
potential Vt, (2) the LG laser beam potential Vlg, and 
(3) the nonlinear mean-field interactions V\. Explicitly, 
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,3 are the trapping frequencies, and fi r(t*) = 
= ^ exp ( — ^ J is the common factor 


of the LG laser potential Vlg? where we have included 
the Stark shift potential, the Raman coupling, and the 
detuning terms. The nonlinear interaction coefficients 
are denoted as g aa f = Airh 2 N a s aa , / M, where N is the 


total number of atoms and a s aa , is the s-wave scattering 
length between the two spins a and a '. 

Similar to our previous work [44], a unitary transfor¬ 
mation Vq/t e^Vq/t gives 
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The second term in the above rotated Hamiltonian is 
an effective SOAM coupling of the form ^ L z cr z , where 
a = H describes the coupling strength and L z = — ihd$ 
is the ^-component OAM operator. One can see explicitly 
the radial dependence oc r~ 2 of SOAM coupling, which 
constitutes a key ingredient in our 2D model. 

The extrema of the LG beam potential along the radial 
direction occurs at R = yjl/2w. Below we use R, k r = 
1/R, and E r = h 2 k 2 /2M as the length, momentum, and 
energy units, respectively. The LG beam potential can 
thus be written in a dimensionless form 


4 lg = ne~ lr \ 21 
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where Q = 4A ° £; is the only dimensionless variable for 
tuning Vlg fixed l. In addition, we consider a har¬ 
monic trap uj x = uj y <C oo z that strongly confines the BEC 
along the z direction, so the system is well described by 
a 2 D GP Hamiltonian. The 2 D trapping potential can 
be written in a dimensionless form as V trap (r) = \Kr 2 , 
where K = M 2 uj 2 1 2 w 4 / 2h 2 is the dimensionless trapping 
strength. In the next section, we will study the single¬ 
particle physics with the dimensionless expression and 
present the ground state diagram as a function of LG 
beam as well as the external trapping potentials. 


III. NON-INTERACTING CASE 


Without interactions, V\ = 0, and the GP Hamilto¬ 
nian becomes a Schrodinger one. Due to cylindrical sym¬ 
metry, the eigenstates can be described by an angular 
quantum number m as well as a radial one n and can be 
written in a form of ( x^ ?n (r)e 2mt<?!) , Xl,n(r)e irn ^ ) T with 
/4 = wz 4 = l [46] . The angular number m is an inte¬ 
ger because OAM quantization of each spin component 
requires m ^ to both be integers. The OAM difference 
— mq = —21 comes from the effect of the SOAM cou¬ 
pling. The radial number n is a non-negative integer due 
to the radial confinement and describes the number of 
radial nodes of the wave function. For a given m, the 
radial part of the wavefunctions T n (r) = ( yq,n Xi,n ) 
and corresponding eigenenergies E n ^ m can be obtained 
by solving a radial Schrodinger equation, 


FT = 


- -d r (rd r ) + |r 2 + | 

T 2 7 


m 2 0 

0 m\ 




~lr\2l 


1 1 
1 1 




(8) 










3 


(a) 


20 
15 

a 10 

5 
0 

60 70 80 90 100 |fl| 


, s " />?'„. 
* / vr 


(b) 


s, 


9\\ 

f/o 

\\\ 

! 



\A'~ 

--A/ 

\ 

f — * / 



k>- < 

Y-O'' 


-3-2-1 0 12 3 



K-\Q\ plane. Different states denoted with OAM quantum 
numbers separate by dashed lines, (b) Energy spec¬ 

trum of the lowest band n = 0 as a function of quasi-OAM 
mn at select points (square, triangle, and circle) in (a). The 
filled symbols denote the ground states, (c) Critical value 
of Raman coupling from double-minima to single-minimum 
transition for blue (A > 0) and red (< 0) detuned lasers and 
for different trapping strengths K. 


where we assume <5 = 0. We find that J5 n?m always in¬ 
creases with n, so it is convenient to interpret the energy 
spectrum by regarding n as a band index. The ground 
state thus lies in the lowest band of n = 0. If we treat 
m as a quasi-OAM number, the lowest band structure 
Eo ?m can exhibit either double minima atm = ±mo with 
0 < mo < l or a single minimum at m = 0, determined 
by O and K. 

In Fig. 2(a), we plot the ground state diagram in \Cl\-K 
plane for the case of / = 2 and A < 0. Each region is 
labeled by the ground-state OAM numbers (mq,mq) for 
each spin species. At a fixed iC, the system undergoes 
transitions from the degenerate ground states m = ±2 to 
m = dzl, and to the single one m = 0 as \Q\ increases. We 
also find that the spin polarization has the same sign as 
m and becomes zero if m = 0. Both characteristics show 
similar physics to that in SLM coupled BEC or SOAM 
coupled ring BEC. Figure 2(b) shows the lowest band 
for selected points at K = 5 in Fig. 2(a). One can see 
how the band structure evolves between different ground 
states. At a fixed £2, we find that K competes with Q and 
moves the system from the single-minimum region toward 
the double-minimum ones. This trend can be understood 
as a result of energetic competition. The external trap 
tends to concentrate the BEC around the central region 
to reduce the potential energy, but rotational energy from 
the centrifugal barrier r~ 2 m ^ in Eq. (8) strongly in¬ 
creases unless the OAM number vanishes. As a result, 
the system favors the majority spin component having 
a low or zero angular momentum and localizing around 
the trap center, even with a certain energetic penalty 
from the minority component and from SOAM coupling. 
In Fig. 2(c), we plot the critical strength £2 C of SOAM 


coupling for the transition between degenerate and non¬ 
degenerate ground states as a function of l at various K 
and different signs of A. We see that fl c increases with 
increasing /, increasing K , or blue-detuned A > 0, con¬ 
sistent with the same physical picture of BEC’s avoiding 
higher rotational energy. Note that the blue-detuned A 
centers the BEC more because the Stark shift creates a 
ring-shaped barrier around the high intensity region of 
the LG beams. 

Before studying the interaction effect, we comment on 
an interesting case, where the difference between two LG 
beams’ azimuthal indices l\ — I 2 is an odd number. In 
this case, m ^ (l\ — h)/2 indicates that m should 

be a half-integer to satisfy OAM quantization. As Q 
increases from zero, the magnitude of m suppresses from 
\(h — Z 2 )/2| until m = ±1/2. Therefore, the lowest band 
always has a double-minima structure and the system 
always has two degenerate ground states with opposite 
spin polarization. 


IV. INTERACTING CASE 
A. Phase diagram 

In this section, we focus on a realistic 87 Rb system, 
compute the full phase diagram, and characterize each 
individual phase as well as various phase transitions. We 
consider two atomic internal states \ ^) = \F = l,rriF = 
0) and | 1) = \F == 1 ,mp = —1) of 87 Rb atoms that 
are coupled through the Raman transitions induced by 
the two LG laser beams with 1 = 2. The radial inten¬ 
sity waist of the LG beam is w = \2\f2fim , which gives 
R ~ 17/ira and E r = 2irh x 0.2 Hz. The ^-direction 
trapping frequency lj z = 2tt x 600Hz is set large enough 
such that the system has an effective 2D geometry [50], 
while the transverse trapping frequency uj x = w y = 1.5 ~ 
6.0 Hz is set isotropic. The s-wave scattering lengths 
are a ^ = 100.86a# and = 100.4ag for 87 Rb, 

where ap is the Bohr radius. In a case of particle num¬ 
ber N = 10 4 , the reduced 2D nonlinear interaction en¬ 
ergy Ei E r is much larger than the single-particle 
rotational energy. Such strong nonlinear interaction will 
drastically modify the single-particle physics, as we will 
see. Note that with this setup, the three nonlinear co¬ 
efficients satisfy g^ga > < 7 ^, which indicates that a ro¬ 
tating spinor BEC without the Raman dressing should 
be miscible [51]. However, we shall see that the wave- 
functions of opposite spins can be spatially separated or 
immiscible in the presence of SOAM coupling. 

We obtain the ground state by solving the GP equa¬ 
tion using the imaginary-time evolution. Figure 3 shows 
the phase diagram in the Raman coupling and detun¬ 
ing Q-S plane, where the colors denote the spin polar¬ 
ization. For comparison, we present the non-interacting 
and interacting phase diagrams in Figs. 3(a) and 3(b), 
respectively. We see that there are four different phases 
in the interacting case according to the spin polarization 
(cr z ) and OAM of the two spins (m^,m^): (i) A stripe 
phase with a density modulation in the angular direc- 
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FIG. 3. (Color online) Phase diagram of the SOAM coupled 
BEC in Q-5 plane for (a) non-interacting system and (b) re¬ 
alistic BEC. The colors defined in the bar graph represent 
spin polarization {cr z }. A stripe phase with small polarization 
occurs at the bottom of (b). The other phases are labled by 
their OAM numbers of the two spin components (mf, 774 ). (c) 
Spin polarization as a function of Raman coupling Q for de¬ 
tuning S — —0.2 ,0.15, and 0.3 E r . (d) Zoom-in of the dashed 
region in (b) showing the stripe phase. Other parameters are 
l o x = uj y = 1.5Hz, N = lx 10 4 , and l — 2. 


tion that is the superposition of (—4, 0) and (0,4) states; 
(ii) a spin-polarized phase, either (—4,0) or (0,+4); (iii) 
a spin-polarized phase, either (—3,+1) or (—1,+3); and 
(iv) a spin-balanced phase (—2,2). We comment that 
the energy difference between these phases can be small 
(compared with typical trapping potential energy, SOAM 
coupling energy, and kinetic energy), and BEC may end 
up with a metastable state in our numerical simulation. 
To obtain the ground-state phase boundaries accurately, 
we start from several initial states to do the imaginary- 
time evolution and choose a final state with the lowest 
energy as the true ground state of the interacting BEC. 
The convergence tolerance of the GP calculation is cho¬ 
sen as lO -4 ^ such that the energy difference between 
the ground state and any metastable state is clearly dis¬ 
cernible. We also emphasize that the ground state of the 
system is robust against small perturbations due to both 
the energy gap from other metastable states and OAM 
quantization. Compared with the single-particle phase 
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FIG. 4. (Color online) The derivative of ground state energy 
dE/dQ as a function of Raman coupling strength £7 for (a) 
S = —0.2E r , (b) S = 0.15 E r , (c) S = 0.3 E r . Compared with 
Fig. 3(c), we see the phase transition is first order where the 
energy derivative is discontinuous whenever the spin polariza¬ 
tion jumps. 


diagram of Fig. 3(a), we see that the interaction shifts the 
phase boundaries and induces a stripe phase [Fig. 3(b)] 
that is described by a state with density modulations. 
Different from the SLM coupling case, the density mod¬ 
ulation here is along the angular direction. In Fig. 3(c), 
we plot the spin polarization as a function of Raman cou¬ 
pling for various detunings. Figure 3(d) is the enlarged 
stripe phase where the spin polarization is zero or small. 
The reason why the stripe phase is not symmetric around 
(5 = 0 is because of the scattering length difference of the 
two spins, which contributes an effective detuning, sim¬ 
ilar to the SOAM coupled ring BEC [44]. Here we see 
the premium detuning value is about S = +0.15EV. In 
Fig. 4, we plot the energy derivative dE/dQ as a func¬ 
tion of Raman coupling. Compared with Fig. 3(c), we 
see that the energy derivative is discontinuous whenever 
the spin polarization jumps or OAM changes, indicating 
that the phase transitions are first order. 

We shall see in the next section that each phase can 
exhibit richer structures in the spin-density distribution, 
including a phase separation in a large parameter region: 
the two spins can be immiscible due to SOAM coupling 
and interaction. In this case, the wave function of the 
two spins is not exactly of the form x cr (r)e 2m<T< ^ ) as in the 
non-interacting case. However, one can still find a set 
of integer winding numbers ( 774 , 774 ) by integrating the 
phase of each spin’s wavefunction over a path enclosing 
the system. Such winding numbers satisfy 774—774 = —21 
and lock with the spin polarization. Therefore, similar 
to the OAM numbers in the non-interacting case, they 
are good quantum numbers for distinguishing different 
phases in the presence of interaction, although the system 
may show different types of phase separation at a given 


B. Individual phases 

In this section, we discuss the individual phases along 
the S = +0.15E r line in the phase diagram of Fig. 3(b) 
in more details. 

For very week Raman coupling and small detuning, 
we see that the system prefers a stripe phase as shown in 
Fig. 5(a). The two spins have density modulations in the 
angular direction, which arise from the fact that there 
are two degenerate dressed states with similar weight 
at quasi-OAM m = — 2 and m = + 2 [green squares in 
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FIG. 5. (Color online) Spin-resolved density distributions of up (upper panel) and down (lower) spins at Raman coupling 
strengths (a) 10F r , (b) 70 E r , (c) 112.Fr, (d) 116F r , (e) 220 E r , and (f) 950F r . The top two rows show the real space density 
distributions, the middle two do the wavefunction phases, and the bottom two do the momentum distributions. The Raman 
detuning is S = +0.15F r . The displayed range of the axis is [— 57/ira, +57/ira] in real space and [— 6k r ,6k r ] in momentum 
space. All other parameters are the same as Fig. 3. 


Fig. 2(b)] that interfere with each other in the real space. 
Because of the superposition of two dressed states, the 
phase [Arg(\h)] of the two spins in the stripe phase is 
almost uniform and their momentum space distributions 
(|<F| 2 ) are centered at k = 0. 

With the increase of the Raman coupling, there is a 
competition between the stripe phase and the occupation 
of a single dressed-state phase. In the SLM coupling case, 
the stripe phase occupies a small parameter region and 
the transition to the single dressed state occurs around 
Q ~ 0.2 E r . However, here we see that the critical value 
for this transition is about H ^ lhE r in the SO AM cou¬ 
pling system. When the particle number or the angular 
momentum l is increased, this parameter region may be 
even larger. Note that the stripe phase has never been 
observed experimentally yet in the SLM coupling system 
because of its small parameter region and the stringent 
requirement to resolve the density modulations separa¬ 


tions (which are at the order of A ^ 1 urn). It might 
be easier to resolve the density modulation separations 
(which are at the order of R ~ 10/im) in the SOAM 
coupled BEC. 

In Fig. 5(b), we show a single dressed-state phase, 
where one of the spin states does not carry any OAM 
and the other spin carries the largest OAM, i.e. m ^ = 0 
and m; = 4. Because of their different OAM, the density 
of up spin is of a rotationally symmetric Gaussian shape, 
while that of down spin is of a ring shape with a hollowed 
center. Similar features are obtained in the momentum 
space density distributions. 

When the Raman coupling is strong enough, the two 
spins may both carry OAM as shown in Figs. 5(c)-5(e), 
where we have (ra^raq) = (—1, +3) or (—2,+2). This 
transition can be understood from the single-particle 
physics picture, where the lowest-energy band evolves 
from a double minima region to a single minimum re- 
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FIG. 6. (Color online) Density distribution of spin up (upper panel) and spin down (lower) for different trapping frequencies 
(a) 1.5Hz, (b) 2Hz, (c) 3Hz, (d) 4Hz, (e) 5Hz, and (f) 6Hz at S — 0 E r and H = 2501+. The displayed range of the axis is 
[—57jum, +57 pm] in real space. All other parameters are the same as in Fig. 3. 



FIG. 7. (Color online) Real-space density distribution and phase of spin up (upper panel) and spin down (lower panel) for LG 
lasers of different OAM (a) l = 2, (b) l = 3, and (c) l = 4 at S = 0 and Q = 2501+. For each figure, the left side is the real-space 
density distribution and the right side is the phase of the wave function. The displayed range of the axis is [—57 jum, +57 jum] 
in real space. All other parameters are the same as in Fig. 3. 


gion. However, an interesting result shows that the two 
spin components appear to be immiscible or spatially 
separated, i.e., wherever one spin dominates, the other 
spin is strongly suppressed. The phase separation oc¬ 
curs in a large parameter region of Q = 110 ~ 8001+ in 
Fig. 3(b) when both spin components carry OAM. The 
emergence of phase separation can be understood from 
the perspective of lowering nonlinear interaction energy 
~ f f(pf + pi) 2 dxdy. When the two spins both carry an¬ 
gular momentum and are not separated, the total density 
of the system must have a hollowed density at the trap 
center. The total nonlinear interaction energy of this sys¬ 
tem is usually larger than the separated phase where the 
total density distribution does not have a hollowed center 
but instead is in a Gaussian shape [51, 52]. 

Among these immiscible states, the two spins carry¬ 
ing opposite OAM +2 and —2 display very interesting 
density distributions, as shown in Figs. 5(d) and 5(e): 
the two holes and two peaks of the two spin components 
could appear in a dipole (d) or quadrupole (e) form. We 
also verify that there are more complicated and rich sep¬ 
arated phases for larger l. The momentum distributions 
of each of the two states are quite different and may be 
distinguished in the time-of-flight experiments. 

In Fig. 5(f), we plot the density distributions of the 
two spins for a very large Raman coupling 12 = 950 E r . 
When the LG beam intensity is too strong, 12 > 8001+, 


both components are trapped around the potential min¬ 
imum of Vlg and forced to mix with each other. Note 
that the Raman coupling and the Stark shift form an 
effective ring-shaped optical trap, which is even deeper 
than the harmonic trap center for such a larger Raman 
coupling strength. In this case, the two components ex¬ 
hibit a miscible ring distribution and their momentum 
distributions are also of a ring-shape, which returns to 
the ring limit [44]. 

V. HARMONIC TRAP AND OAM INDUCED 
PHASE TRANSITION 

In this section, we study the effects of the harmonic 
trapping potential and larger OAM number of the LG 
beams. From the single-particle results in Fig. 2(c), we 
see that the trapping strength may drive a quantum 
phase transition between the states with different OAM. 
In Fig. 6, we plot the real-space density distributions 
of the two spins with increasing trapping potential from 
1.5Hz to 6Hz. The two spins are pushed inward when the 
trapping potential is increased and the OAM of the two 
spins changes from (—2,2) to (—3,1) and then (—4,0), 
which is consistent with the energetic competition de¬ 
scribed in Sec. III. In previous SLM coupled systems, 
the trapping strength was relatively weak and thus there 
were no observable effects with the change of trapping 
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frequencies. Our results demonstrate that the SOAM 
coupled BEC has a significant response to the trapping 
potential. 

In Fig. 7, we plot the density distributions and phases 
of the two spins for (a) Z = 2, (b) Z = 3, and (c) l = 4. We 
see that the two spins are in a (—2, 2) separated phase for 
l = 2, then in the (—6,0) spin-polarized miscible phases 
for l = 3, and eventually a stripe phase from the superpo¬ 
sition of (—4, 0) and (0, +4) dressed states for l = 4. This 
can be understood from the l dependence of the critical 
Raman coupling strength, similar to the single-particle 
results in Fig. 2(c). 

VI. CONCLUSION 

In conclusion, we studied the quantum phases of a 
SOAM coupled BEC with realistic experimental parame¬ 
ters in a 2D geometry. The phase diagram shows distinct 
physics from a SLM coupled BEC due to the naturally 
strong interaction and OAM quantization, and even from 
a SOAM coupled ring BEC due to the radial inhomogene¬ 
ity. The interaction-induced stripe phase and different 


immiscible states have been characterized by the spin- 
resolved position and momentum distributions. Various 
quantum phase transitions induced by the Raman cou¬ 
pling, detuning, trapping potential, or LG beams’ OAM 
have been reported too. Although the results are pre¬ 
sented for 87 Rb atoms, our model and treatment can be 
easily applied on others such as 23 Na. For a typical Ra¬ 
man coupling strength of 1 kHz and an experimental time 
scale of 1 s, the heating effects in both 87 Rb and 23 Na 
systems should be weak [28, 44]. 

Since the dynamics of ultracold atoms under the ap¬ 
plication of LG beams have been experimentally inves¬ 
tigated by several groups [53-58], and ongoing experi¬ 
mental effort is along this direction for the realization of 
SOAM coupling [59], the achievement of a BEC ground 
state with SOAM coupling would be expectable in the 
near future. Our results hence provide timely predictions 
for experimental observation. 
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